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Recent years have seen the discovery of a chiral spin liquid state - a bosonic analogue of a fractional Quantum 
Hall state first put forward by Kalmeyer and Laughlin in 1987 - in several deformations of the Heisenberg model 
on the Kagome lattice. Here, we apply state-of-the-art numerical techniques to one such model, where breaking 
of the time-reversal symmetry drives the system into the chiral phase. Our methods allow us to obtain explicit 
matrix-product state representations of the low-lying excitations of the chiral spin liquid state, including the 
topologically non-trivial semionic excitation. We characterize these excitations and study their energetics as the 
model is tuned towards a topological phase transition. 


I. INTRODUCTION 

For decades, researchers have sought to find examples of 
spin liquid phases [1], i.e. Mott insulators that are at low¬ 
est temperatures not characterized by local order parameters. 
While numerical simulations are frequently invoked to ex¬ 
plore these phases, the reliable characterization of spin liq¬ 
uid phases has remained a great challenge. Recent years have 
seen a fast evolution of numerical methods that have facili¬ 
tated a paradigm shift away from characterizing spin liquid 
phases by the lack of local order towards a characterization 
through positive, universal signatures. This is particularly the 
case for gapped, topologically ordered spin liquids, where an 
extensive theoretical framework exists for the universal prop¬ 
erties of these phases [2]. 

A recent proposal to numerically extract the behavior un¬ 
der modular transformations [3] has, where practically feasi¬ 
ble, allowed an unambiguous identification of the universal 
properties of a topological phase. In Ref. 4, it was demon¬ 
strated that this technique can be applied to matrix-product 
state approximations [5-7] of the ground states of some sys¬ 
tem on thin, infinitely long cylinders. This allows the extrac¬ 
tion of both the entanglement spectrum [8] and the modular T 
and S matrices that fully characterize topological bulk excita¬ 
tions using unbiased numerical techniques for a large class of 
Hamiltonians. 

Here we extend these methods and use them to construct ex¬ 
plicit wavefunctions for isolated topologically non-trivial and 
trivial excitations. We characterize these excitations through 
their SU(2) and momentum quantum numbers and explore 
their energetics, allowing us to gain understanding of the ex¬ 
citation spectrum of the topological phase. We apply this 
method to an SU(2)-invariant spin-1/2 model on the Kagome 
lattice, where a recent study [9] using the aforementioned 
techniques established the presence of a topologically or¬ 
dered chiral spin liquid phase [10]. This phase emerges when 
time-reversal symmetry is broken in the well-known Kagome 
Heisenberg antiferromagnet by adding a sufficiently strong 
coupling to the scalar spin chirality [11, 12]. 

Beyond being able to gain access to explicit wavefunctions 
for topological excitations in a spin system, our methods al¬ 
low us to study the interplay of topological order and SU(2) 
symmetry in this particular realization of the chiral spin liq¬ 
uid. In a model that both has a microscopic symmetry and 


emergent topological order, a possible scenario is that of sym¬ 
metry fractionalization, where the anyonic excitations carry 
fractional charge of the microscopic symmetry, as opposed to 
integral charge carried by conventional local excitations. Fa¬ 
mously, this is the case in the quantum Hall effect [13], where 
the anyons carry fractional electric charge. It is thus natural to 
expect that the excitations of the chiral spin liquid, whose uni¬ 
versal properties are described by a bosonic iz = 1/2 Laughlin 
state [10], carry fractional SU(2) charge [14]. We will demon¬ 
strate that this is indeed the case in this particular model. 

Given our numerical characterization of the low-lying ex¬ 
citations, we are able to establish a close connection between 
topological and non-topological degrees of freedom deep in 
the chiral phase, where we show that the lowest spin excita¬ 
tions can be identified as pairs of non-interacting fractional 
excitations. Tuning the model towards a topological phase 
transition out of the chiral phase, we observe the onset of in¬ 
teractions between the fractional excitations. We furthermore 
observe a major reorganization of the low-lying singlet exci¬ 
tations which appears to drive the transition out of the chiral 
phase. 

The paper is structured as follows: In Section II, we intro¬ 
duce the model that will be studied. In Section III, we describe 
an ansatz that allow us to obtain explicit wavefunction rep¬ 
resentations for excitations for infinite quasi-one-dimensional 
systems. In Section IV, we discuss our numerical results. 


II. MODEL 

Following a number of exact parent Hamiltonians [15-18] 
as well as realizations in topological flat-band models [19- 
21], recent years have seen the chiral spin liquid realized in 
a number of deformations of the Heisenberg model on the 
Kagome lattice [9, 22, 23]. The model we focus on here is 
given by [9] 

H = cos{9) ^ Si ■ Sj + sin(0) ^ Si ■ {Sj x Sk), (1) 
fj) 

where S denotes spin-1/2 operators acting on the sites of a 
Kagome lattice, and the sum in the first term runs over all pairs 
of nearest neighbors, whereas the second sum runs over the el¬ 
ementary triangles of the Kagome lattice and sites are ordered 
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clockwise around each triangle. The first term describes the 
well-known Heisenberg antiferromagnet, which has received 
much attention over the years [24-26]. Recent evidence from 
large-scale density-matrix renormalization group calculations 
has favored a gapped phase without local order [27] which ap¬ 
pears topologically ordered with a topological entanglement 
entropy of log(2) [28, 29]; at the same time, calculations using 
variational Monte Carlo hint at a gapless state [30, 31]. As¬ 
suming the state is indeed gapped and topologically ordered, 
the value of the topological entanglement entropy would be 
consistent with two distinct topological phases, which we will 
refer to as ’’toric code” [32] and ’’double semion” [33] phases, 
although their history as distinct Z 2 gauge theories predates 
the toric code. While these phases cannot be distinguished by 
their edge or their ground state degeneracy, their bulk exci¬ 
tations exhibit different exchange statistics. Which of these 
phases is realized in a putative time-reversal symmetric topo¬ 
logical phase in the Kagome Heisenberg antiferromagnet is at 
this point an unresolved question; recent work has put forward 
energetically favorable ansatz states for both candidates [34- 
37], but has also placed severe symmetry constraints on the 
double semion wavefunction [38]. 

The second term of Eqn. (1), coupling to the scalar spin 
chirality [ 11 , 12 ], breaks time-reversal while preserving the 
SU(2) spin symmetry and leads to the formation of the chiral 
phase for 9 > O.OStt [9]. The universal properties of this chi¬ 
ral phase are captured by a bosonic z/ = 1/2 Laughlin state, 
which exhibits a two-fold degeneracy of its ground state on the 
torus or infinite cylinders, where the states can be character¬ 
ized by different topological flux through the torus or cylinder; 
bulk excitations that obey semionic statistics, i.e. the wave- 
function acquires a phase factor i when two quasiparticles 
are interchanged; and a chiral gapless edge state described by 
an SU(2)i Wess-Zumino-Witten (WZW) model with a chiral 
central charge of c = 1. This gapless edge is observable not 
only in the energy spectrum when placed on a system with a 
boundary, but also in the entanglement spectrum [ 8 ]. Other 
aspects of this phase have recently been studied using a vari¬ 
ety of numerical [39] and analytical [40, 41] techniques. 

No signs of time-reversal symmetry breaking have been ob¬ 
served for the Heisenberg model itself. Therefore, we expect 
a transition to take place between the chiral spin liquid phase 
and a putative time-reversal symmetric phase surrounding the 
Heisenberg point. Such a transition has been analyzed re¬ 
cently in a different model in Ref. 42, as well theoretically 
in Ref. 43. For the model considered here, we have shown [9] 
that the transition must take place in the range Q <9 < O.OStt. 
Given the inherent limitations in terms of system size, which 
are accentuated by the increasing correlation length as the 
phase transition is approached, it is challenging to locate the 
transition more precisely and argue whether its nature is first- 
order or continuous. 

The symmetry-enriched properties of this realization of the 
chiral spin liquid can be determined from the entanglement 
spectrum reported in Ref. 9, where we found the degeneracies 
of the entanglement spectrum to agree with the descendants of 
the two primary fields of the SU(2)i WZW model. Within this 
model, the states can be organized into SU(2) multipletts [44]. 


In the tower of descendants of the identity field, the spins for 
the first four levels are 0 , 1,0 0 1,0 0 1 0 1 ; we find this spec¬ 
trum to be realized on the edge of a cylinder with trivial flux. 
Conversely, a cylinder with semionic flux realizes the descen¬ 
dant spectrum of the spin- 1/2 field, where we have 1 / 2 , 1 / 2 , 
1/203/2, 1/201/203/2. We note that these SU(2) quantum 
numbers are emergent and the same level counting is found in 
models that lack a microscopic SU(2) symmetry [9, 19-21]. 
However, for models that do have a microscpic SU(2) symme¬ 
try, we can assign to each level in the entanglement spectrum 
a micoscopic SU(2) quantum number. Performing this iden¬ 
tification for the ground states of ( 1 ), we find that the micro¬ 
scopic quantum numbers coincide with the emergent quantum 
numbers in the CFT spectrum, which implies that the semion 
carries a fractional spin of S' = 1/2. We emphasize that no 
local, topologically trivial excitation can carry such spin; a lo¬ 
cal spin flip, for example, corresponds to an S = 1 excitation. 
Considering the fusion rule of the semion, s 0 s = 1, and com¬ 
paring this with the well-known fusion rules for SU(2) spins, 
e.g. 5 0 I = 0 0 1 , we can generally identify half-integer 
spin excitations as topologically non-trivial and integer spin 
excitations as trivial. 


III. AN ANSATZ FOR LOCALIZED EXCITED STATES 

To access the properties of the model (1) computationally, 
we consider the Hamiltonian on a cylinder with infinite length 
and finite width. A key advantage of the geometry of infinite 
cylinders is that for a topologically ordered system, we can 
obtain the full set of degenerate ground states in a basis of 
states with well-defined topological flux [4]. 

For the Hamiltonian (1) on an infinite cylinder, restricted 
to the parameter range O.OStt < 9 < O.Stt, Ref. 9 finds 
two ground states: |<l>]i) and |$s). Here, |$]l) denotes a 
ground state with identity flux and |$s) has a semion thread¬ 
ing through the cylinder. The Schmidt decomposition of |$i) 
and I $ s) with respect to partitioning the system into two semi¬ 
infinite cylinders reads 

|<i>i)=^A^|Li)|i?i), (2) 

OL 

\<^.)=Y.^l\Ll)\Rl) , (3) 

OL 

where |Lq) (|i?a)) form an orthonormal basis of vectors of 
the left (right) partition of a cylinder. Aq, are the Schmidt co¬ 
efficients and (Atj)^ are the eigenvalues of the reduced density 
matrix of the left (right) semi-infinite cylinder. Fig. 1 gives a 
graphical representation of the Schmidt decomposition of the 
ground states |<i)i) and |<i>s). 

Our ansatz for the excited states is based on the idea of 
replacing the Schmidt coefficients A in (2,3) by a matrix G 
allowing for off-diagonal terms |La)|i?^) to contribute. The 
ansatz for the excitation takes the form 

|'I^)=5lGa/3|Tc.)|i?/3), (4) 

a,/3 
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FIG. 1. (Color online) The Schmidt decomposition of the ground 
states of the Hamiltonian (1) on an infinite cylinder. |<l?i) carries no 
topological flux while |il?s) has a semion s propagating through the 
cylinder. 


where the matrix G is normalized such that 
(tl'I'P) = Tr(G^G) = 1 and is energy optimized un¬ 
der constraint ($|'P) = ^aGaa = 0 to ensure that 

the excitation Ith) is orthogonal to the ground state |$). 
The ansatz (4) represents an excitation that is localized in 
the longitudinal (i) direction of a cylinder and completely 
delocalized in the transversal (y) direction. Specifically, the 
excitation is localized within the correlation length in the x 
direction and has a has a well-defined transverse momentum 
ky. Throughout the remainder of this section, we will use 
the notation to denote different types of excitations, 

where a stands for the topological charge of an excitation, 
which can take the values o = 1 (identity) and s (semion) for 
the case considered here, and y, labels its quantum numbers 
under the local symmetry. 

For an SU(2) invariant ground state, an index a that labels 
the Schmidt vectors \La), \Ra) in (2,3) carries SU(2) quan¬ 
tum numbers; the total spin J and its z component Jz, where 
J = Ri is the total spin operator. This feature allows us to 
access excitations with well-defined SU(2) quantum numbers 
by enforcing further constraints on the matrix G. We will de¬ 
note such excitations as with y = S (singlet, J = 0), 

T (triplet, J = 1) or Z? (doublet J = 1/2). 

For example, a singlet excitation above the ground state 
with the identity flux is obtained by replacing the Schmidt co¬ 
efficients in ( 2 ) by a matrix G^’'^, 

ivi/i.^) = ^G;i’;iLi)ii?i), (5) 

a,/3 

1 S 

where the only non-zero elements G^’^ are such that the in¬ 
dices a and /? fuse into a singlet. Matrix G^’® in (5) is the 
subject of an energy optimization. 

Let us now replace the Schmidt coefficients by a matrix 
G^’^, which is constrained in a way that ^ 0 only if 

indices a and /3 are such that forms a state with 

J = 1 and = 1. The above constraint leads to an ansatz 
for a spin -1 excitation 

=Y.GII \Ll)\Rl) . ( 6 ) 

In the following, will be referred to as triplet excita¬ 
tion. Strictly speaking, is one of the three states of 


spin 0 N 
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(b) 
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FIG. 2. (Color online) An intuitive picture illustrating the ansatz 
Pqj- (Jifferent types of excitations, (a) The Schmidt coeffi¬ 
cients A’^ from Fig. i(a) are replaced by a matrix to form a sin¬ 
glet excitation (b) An ansatz for the triplet excitation 

is obtained by inserting matrix Matrices G^’^ and are 

energy optimized under different set of constraints discussed in the 
text, (c) In turn, an ansatz for the topologically non-trivial excita¬ 
tion is constructed by gluing two semi-infinite cylinders, one 

with identity flux and the other with semion flux inside. The semion 
flux propagating in the right part of the cylinder creates a spin-1/2 
semionic excitation. All the excitations are localized in the shaded 
region and have a well-defined transverse momentum ky. 


an SU(2) triplet. States with Jz = 0 and Jz = —1 may be 
constructed in a similar way. Correspondingly, matrix G^’^ is 
obtained as a result of minimizing the energy of 14/ 

Eqs. (5,6) represent topologically trivial (neutral) excita¬ 
tions above the ground state with identity flux. Having access 
to the ground state with the semion propagating inside the 
cylinder, we can generalize the ansatz (4) to describe topo¬ 
logically non-trivial excitations [45-47]. The ansatz for such 
an excitation takes the form 

(7) 

a,/3 


Here, we force the state to have well-defined semion 

flux inside the right semi-infinite cylinder and identity flux in 
the left part of the cylinder. Semion flux is seen to create a 
spin- 1/2 topological defect around the position of the matrix 
G®’^. In this setting, the matrix G®’^ is chosen to minimize 
the energy, subject to a normalization constraint only. 

Fig. 2 provides an intuitive illustration of the ansatz for the 
different types of excitations described above. Appendix A 
discusses the specific realization of that ansatz in which the 
Schmidt vectors |Lq,) and \Ra) are given by matrix-product 
states. In that setting, we explain how to compute the energy 
and momentum of topologically trivial and non-trivial excita¬ 
tions. Appendix A also describes the computation of overlaps 
between ground states and excited states on an infinite cylin¬ 
der. 
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FIG. 3. Section of an infinite XC8-4 cylinder. Wavy lines denote 
boundary conditions in the y direction of the cylinder. 



Q I'K 

FIG. 4. (Color online) Energy of a single, isolated semionic excita¬ 
tion above the ground state of (1). Results are obtained for a cylinder 
of width 8 with a bond dimension of x = 4096. 


IV. RESULTS 

We have applied the methods discussed in Section III 
- specified in the context of matrix-product states in Ap¬ 
pendix A - to the model given in Eqn. (1). Throughout this 
work, we use a cylinder in the XC8-4 geometry, as shown in 
Fig. 3 and quote results obtained for a bond dimension of up 
to X = 4096. 

In Fig. 4, we show the energy of a single semion, 
throughout the phase diagram of Hamiltonian (1). We observe 
that the semion energy increases as the model is tuned from 
the fully chiral point at 0 = 7r/2 to the Heisenberg point 0 = 
0. For all values of 0, the semion energy decreases slowly 
as the bond dimension is increased (in the figure, note that 
the bond dimensions are spaced exponentially). Crucially, the 
semion energy seems to remain finite throughout the entire 
phase diagram, including at very small values of 0 in the range 
0 < 0 < O.OStt where we expect a phase transition. 

In Fig. 5, we show, in addition to the energy cost of two 



Q j'K 

FIG. 5. (Color online) Energies of low-lying topological and trivial 
excitations on a XC8-4 cylinder with x = 4096. Excitation energies 
shown here are measured in the identity sector; results for the semion 
sector are very similar. 

semions at infinite distance, 2i?®®™, the variational bounds on 
the energies of local, topologically trivial excitations of spin 
J = 0 (singlet) and J = 1 (triplet), which we will below 
denote as A^® and ®, respectively. Fach of these quantities 
can be extracted for excitations above the ground state in the 
identity and semion sectors; our results show that these gaps 
agree within reasonable numerical accuracy and we will thus 
focus on the gaps above the identity sector in the following. 

A key observation is that at the fully chiral point 0 = 7r/2, 
where the first term of Fqn. (1) - the Heisenberg coupling - 
vanishes, the energies of a local singlet and triplet excitation 
approximately coincide with twice the energy of a semion, 
that is 

A^^® « a!)?® « 2 F;®®“. (8) 

This suggests the interpretation that at least for 0 = 7r/2, the 
semions interact very weakly at short distances and the local 
singlet and triplet excitations, which are topologically trivial, 
consist of a pair of semions. This is consistent with the two 
possible SU(2) fusion outcomes for a semion of spin J = 1/2, 
given by the usual SU(2) spin addition rule 1/201/2 = 0©1. 
In the following, we will examine this assertion in more detail 
and in particular carefully examine whether it extends away 
from 0 = 7r/2. 

In the data shown in Fig. 5, we see that as the Hamiltonian is 
tuned away from the purely chiral point, the energies of local 
excitations behave markedly different. Focusing first on the 
behavior of the energies for topologically trivial excitations 
A^® and A!^®, we find that they remain finite throughout the 
phase diagram. This is to be expected even in the case of 
a continuous phase transition from the chiral spin liquid to 
a topologically distinct gapped phase around the Heisenberg 
point as the transition is rounded to a crossover due to the 
finite size of the system, and the extent of the excitation being 
considered is limited due to the finite bond dimension of the 
MPS. We also note that the upper bounds on the excitation 
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FIG. 6. (Color online) Interaction energies extracted using the 
ansatz (9). Results are obtained using x = 4096. 


gaps given by and shown in Fig. 5 do not seem to be 
tight, as we can construct a (non-local) excitation of two well- 
separated semions, which, assuming that interactions between 
anyons decay with distance, will have energy 

However, the data also clearly shows that while the triplet 
gap Ay'^ remains quantitatively unchanged across the transi¬ 
tion, the singlet gap A'j"^ drops signihcantly as the transition 
is approached. This would be consistent with a closing of 
the singlet gap in the thermodynamic limit, while the triplet 
gap remains hnite across the transition even in the thermody¬ 
namic limit. Unfortunately, reliable extrapolation to the ther¬ 
modynamic limit in critical region is beyond the capabilities 
of current DMRG methods. 

The interpretation of the local singlet and triplet excitation 
as a pair of semions naturally suggests an interpretation of the 
difference between the energy of the trivial excitation and that 
of two semions as an interaction energy between semions. 


TT-int 

ZT'int 

H/rp 



- 


(9) 


This can be thought of as the energy cost of bringing two in- 
hnitely distant semions close to each other (on the scale of 


FIG. 7. (Color online) Singlet (top panel) and triplet (bottom panel) 
excitation energies for x = 4096. Here, the symbols denote the ky 
momentum of the excitation, and the thickness of the lines indicates 
the overlap between excited states. 


the correlation length), or conversely as the energy gain from 
locally creating a pair of semions and bringing them to a dis¬ 
tance from each other that is much larger than the correlation 
length. When i?™* is positive, the two semions thus interact 
repulsively, while they interact attractively for < 0. 

Our numerical results for this quantity are shown in the top 
panel of Fig. 6. At 0 = 7r/2, we have that i?™* « 0 both 
in the singlet and the triplet channel. For 7r/4 < 9 < 7r/2, 
we hnd that the interaction energy becomes repulsive, > 
0, in both channels. As the transition regime is approached 
further, we observe drastically different behavior in the singlet 
and triplet channels: while the interaction energy in the triplet 
channel remains qualitatively unchanged as the transition is 
approached, the interaction energy in the singlet channel drops 
and eventually changes its sign at 0 « O.IStt. Note that in 
Fig. 5 this behavior is reflected in the fact that Ag'^ crosses 
2 ^sem around that value of 9. 

The sudden crossover in suggests a further investiga¬ 
tion of the question whether the interpretation of the localized 
singlet excitation as a pair of semions remains valid around 
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0 « 7r/4. To this end, we have calculated the hdelities of 
low-lying singlet and triplet excitations, as shown in Fig. 7. In 
the hgure, data points indicate the energy of low-lying local 
excitations, and the thickness and shading of the lines con¬ 
necting them indicate the hdelity of a pair of excitations. In 
the triplet sector, we find that the lowest excitation, which car¬ 
ries transverse momentum ky = 0 forms a branch of excita¬ 
tions with high overlap that is well-separated from the rest 
of the spectrum, both in terms of its energy as well as its h- 
delity. This strongly suggests that the interpretation of the 
triplet as a pair of semions holds not only at the fully chiral 
point 6 — TT12, but rather holds throughout the entire phase 
diagram up to the transition. In contrast to this, the spectrum 
of excitations in the singlet sector shows a major reorganiza¬ 
tion around 9 = 7r/4, coincidental with the regime where the 
excitation energy changes drastically. In particular, the branch 
of the lowest excitation for 0 > 7r/4 moves up in the spectrum 
for 0 < 7r/4. This must be taken as indication that the inter¬ 
pretation of the local singlet excitation as a pair of semions 
must become unreliable in the vicinity of 6* = 7r/4. 

We can further resolve this data by calculating the vari¬ 
ational energies of the ansatz states independently for the 
Heisenberg and chiral interaction, which we will now denote 
as and respectively. These quantities are shown in 
the middle and bottom panels of Fig. 6. The results show that 
the interaction energy due to the chiral term for the triplet ex¬ 
citation is very small throughout the phase diagram; in the sin¬ 
glet channel, the interaction energy increases sharply around 
6 — 7r/4. For the interaction due to the Heisenberg term, 
we observe that it remains qualitatively the same in the triplet 
case throughout the entire phase diagram, while in the singlet 
case it shows a sharp decrease around 6 = ir/A, where the 
Heisenberg term becomes the dominant term in the Hamilto¬ 
nian. This is also consistent with the major reorganization of 
the excitation spectrum at 9 = 7r/4 driven by the presence of 
a Heisenberg of comparable strength to the chiral term. 


V. CONCLUSION 

We have described matrix-product state techniques to study 
the excitations of a topological phase on inhnite cylinders. 
These techniques allow us to obtain both topologically triv¬ 
ial excitations as well as non-trivial excitations, which are ob¬ 
tained as domain walls between cylinders of different topolog¬ 
ical flux. In our approach, the local excitations are confined to 
within roughly one correlation length of the MPS, and ener¬ 
gies obtained in this way thus provide an upper bound on the 
energy of the excitations. These bounds could be improved by 
allowing larger spatial extent of the excitations. 

We apply these techniques to a Heisenberg model where 
a chiral three-spin term breaks time-reversal symmetry ex¬ 
plicitly and favors a chiral spin liquid phase. The model we 
study interpolates between the chiral spin liquid in the limit 
of sufficiently strong time-reversal symmetry breaking, and 
the time-reversal symmetric phase that surrounds the unper¬ 
turbed Heisenberg point. We show that in the limit where 
only the three-spin chiral term is present, the low-lying sin¬ 


glet and triplet excitations are given as pairs of weakly inter¬ 
acting semions. As the system is driven towards the Heisen¬ 
berg point, the semions begin to interact and the energetics 
of the low-lying excitations change significantly. We observe 
that for the finite systems available to our numerical methods, 
both the energy of a single semion as well as the energy of 
a local triplet excitation remain finite across the entire phase 
diagram. The energy of a local singlet excitation, on the other 
hand, drops sharply as the transition is approached, which is 
accompanied by a major reorganization of the low-lying sin¬ 
glet excitations in the chiral phase. 
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Appendix A: Matrix-product ansatz for localized excited states 

In this Appendix we discuss a matrix-product state ansatz 
for localized excited states outlined in Section III. 

The system under consideration consists of a cylinder with 
infinite length and finite width W, as measured by the lat¬ 
tice unit cells. Fig. 3 in Section III shows the Kagome lattice 
with our choice of boundary conditions in the transverse di¬ 
rection for W = 4. In this setup, the Hamiltonian remains 
translationally invariant in both longitudinal and transverse di¬ 
rections. However, a matrix-product ansatz for the state on 
an infinite cylinder is translationally invariant in x direction 
only. Translational invariance in the y direction will only be 
approximately recovered through energy minimization. The 
MPS ansatz considered here for ground states jfl)) and excited 
states |T') consists of a unit cell of n different tensors that is 
repeated throughout an infinite cylinder. We denote the MPS 
tensors as {Ai} and {Bi}, where i = 1,... ,n enumerates the 
tensors on inequivalent sites within the unit cell, and denote 
the bond dimension of each tensor as %. Tensors {Ai} ({Bi}) 
correspond to the left (right) semi-infinite part of a cylinder, 
see Fig. 8. 

The Schmidt decomposition of the ground state 1$) on an 
infinite cylinder takes the form 

1$) , (AI) 

a. 

where the Schmidt vectors |Lq,) and li?^) are represented by 
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FIG. 8. (Color online) Matrix-product state representation of a state 
on an infinite XC8-4 cylinder with finite width W = 4 lattice unit 
cells. Green line represents the mapping of an infinite cylinder to 
a Id path. The MPS ansatz is composed of a unit cell of n = 12 
tensors. Boudary conditions are marked with a wavy line 


the MPS tensors (compare with Fig. 8): 


AAA 


\L^) = 

(A2) 


(A3) 

We emphasize that, in (A2), A represents not just 
tensor but the entire unit cell of tensors Af. 

one single 

A Ai A 2 An 

(A4) 

and similarly for tensor B in (A3). In the canonical MPS rep¬ 
resentation of |Lq,) and tensors A and B are related by 

B = A“^AA, where A is the diagonal matrix of the Schmidt 
coefficients in (Al). 

As described in Section III, an excited state is constructed 
by replacing the Schmidt coefficients A in (Al) by a matrix G: 

A A G B B 

IT-) = 

, (A5) 


where G is a x x X bond matrix without physical index. Nor¬ 
malization of |tlf) requires that = 1. The 

matrix G is constrained by ^aGaa = 0 to as¬ 

sure that the excited state |'I') is orthogonal to the ground state 
1$). Forcing additional constraints on the matrix G discussed 
below, leads to an ansatz for a singlet or a triplet 

excitation with corresponding matrices and 

Additionally, a topologically non-trivial excitation is ob¬ 
tained by connecting two semi-infinite cylinders carrying 
identity and semion flux. Following Eq. (7), the MPS ansatz 


is given by 

a1 a1 ^S.D TD S T3 S 

(A 6 ) 

where A^ forms the Schmidt vectors |L^) of the ground state 
having identity flux threading through the cylinder. Sim¬ 
ilarly, cotTesponds to the Schmidt vectors |i?®) of |$s) 
carrying the semion flux. 

In order to build the above localized excitations with SU(2) 
charge J and = J it is enough to restrict our attention 
to the U(l) subgroup generated by J^. Correspondingly, we 
have used an MPS that is made of U(l) symmetric tensors. In 
that case, an index a labeling Schmidt vectors |Lq,) and \Ra) 
in (A2,A3) is represented by a pair a — (c, i), where c is a 
U(l) charge and i labels different Schmidt vectors carrying 
the same charge c. The matrix G in (A5) is U(l) symmet¬ 
ric, which means that it maps charge c to (possibly different) 
charge c + Ac. Technically, this implies that G has a block 
structure in which the only non-zero elements G(^c,i),(c',i') ^6 
such that c' = c — Ac, where Ac is the same for all possible 
values of c. In this setting, the excited state (A5) is seen to 
carry U(l) charge Ac. 

A singlet excitation A) obtained by choosing a block 
structure with Ac = 0 and ensuring that the resulting ex¬ 
cited state does not belong to a larger multiplet. This require¬ 
ment can be verified numerically by computing the expecta¬ 
tion value of (J + Jy + J^, where J is the total 

spin operator. In practice, we only need to evaluate (J on 
a finite, cylindrical region surrounding the matrix G^’^, de¬ 
noted here by A. The result is then compared against (Ja)^ 
computed for the ground state which is known to be a singlet. 
The difference (T'|(Xt)^|T') — gives J(J + 1), 

where J is the total spin of an excitation. The reason to con¬ 
sider the difference instead of (T'|( J^)^|T') is that on a finite 
region A, ('F|(J^)^|'F) contains contributions coming from 
the boundary of A. However, it follows from the construction 
of the excited state (A5), that the same boundary contribution 
is present in (<i)|(Xt)^|$). 

In order to access a = 1 state of a triplet excitation 
matrix is constrained with Ac = 1; then one 
must verify that such an excitation does not belong to some 
larger multiplet, which can be achieved in a fashion analogous 
to the above prescription for the singlet excitations. 

In the construction of a topologically non-trivial excitation 
matrix G®’^ has a peculiar block form. In this case, 
the indices a and /3 of the Schmidt vectors |L^) and |i?|) 
carry SU(2) quantum numbers that fuse to half-integer repre¬ 
sentations of SU(2). In the identity topological sector, |L^) 
carries integer J, while in the semion sector, |i?|) is associ¬ 
ated with half-integer J. Thus, the matrix G®’^ maps U(l) 
charge c to U(l) charge c + Ac, where Ac is half-integer. 
We note that this immediately implies that the semion carries 
fractional U(l) charge (or equivalently, half-integer J). 

To obtain the energy of a local excitation, we measure all 
Hamiltonian terms in a region A that contains the excitation 
and is sufficiently large compared to the correlation length. 
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see Fig. 2(a)-(c). Performing this for the same region both 
in the ground state |$) as well as the excited state j^*), the 
energy of the local excitation is given as the energy differ¬ 
ence. Assuming that the ground state energy was obtained 
accurately, the energy difference introduced above provides 
an upper bound for the corresponding excitation gap. This as¬ 
sumption is reasonable, since DMRG is known to accurately 
capture the properties of the ground state (which may not be 
true for the excited state). In the case of a topological de¬ 
fect, the additional subtlety arises that the two ground states 
|$i) and |$s) in general have slightly different energy density 
(due to a splitting that vanishes exponentially in the circumfer¬ 
ence of the cylinder W). In this case, the energy is compared 
against the average of the energies of the two ground states 
in region A. Note also that in this approach, we can evaluate 
the energy separately for different terms of the Hamiltonian 
(1). Section IV discusses independent contributions due to 
Heisenberg and scalar spin chirality terms. 

In addition, we can obtain the transverse momentum ky, 
i.e. the momentum related to translations around the cylinder, 
of a local excitation. This computation is crucially based on 
the observation that translation Ty by one lattice unit cell in y 
direction is a product of operators Tq acting within one MPS 
unit cell. As can be seen from Fig. 8, the overlap between 
MPS tensors A, A* and the translation Tq is given by 



(A7) 


This observation allows us to measure the expectation value 
of translation Ty, which for a momentum eigenstate should 
equal . The transverse momentum is thus given by 


A A G B B G 



A* A* G* B* B* G* 


(A8) 

As the boundary objects (denoted as P and Q in (A8)) are 
only defined up to a phase, this expectation value is also only 
defined up to a phase. However, since the boundary objects 
are the same for the ground state |<i)) and excited states |T') 
that differs only by a bond matrix G, we can obtain the mo¬ 
mentum of an excited state relative to the ground state, which 
is assumed to have momentum ky = 0. For a well-converged 
excited state, the momentum will be close to one of a set of 
discrete eigenvalues ky = “^m, m = Q,... ,W — 1. If the 
bond dimension x is insufficiently large, we often find that ky 
is not near one of these values, indicating that the state |T') is 
not converged to an eigenstate of the translation operator Ty. 


Finally, let us consider two infinite MPSs 

A A G B B 

10(G)) = ■■■ ■■■, (A9) 

A^ G' B^ 

|0'(G')) = ■■■ (AlO) 

States |0(G)) and |0'(G')) may represent localized excita¬ 
tions or ground states. In the latter case. Gay = Say^a and 
G(j,^ = Say^a’ where Aq, A), are the Schmidt coefficients 
(Al). 

The overlap per unit cell between ground states |0(A)) and 
|0'(A')) is defined as the absolute value g of the dominant 
eigenvalue of the mixed transfer matrix T, 


A 



A'* 


For normalized |0(A)) and |0'(A')), we have g S [0,1]. In 
the case of p = 1, tensors A and A' describe the same state 
and are related by a unitary gauge transformation given by 
the dominant eigenvector of T. Note that the corresponding 
transfer matrix built by tensors B and B' has the same spec¬ 
trum, since tensors A and B as well as A' and B' are related 
by B = A-^AA and B' = (A')”^A'A'. 

For the purpose of this paper we propose a measure of the 
overlap between localized excitations of the Hamiltonian ( 1 ) 
for nearby values of the parameter 9. In this approach we 
assess the overlap locally, in the region surrounding the exci¬ 
tations, and neglect the fact that far away from the excitation, 
states |0(G)) and |0'(G')) have overlap per unit cell strictly 
smaller than 1. 

To this end, we assume that the ground states |0(A)) and 
|0'(A')) have the overlap per unit cell close to 1, which means 
that tensors A and A' are (up to a gauge) almost identical. 
This is the case of the ground states of the Hamiltonian (1) 
for nearby values of the parameter 9, as found in [9] and fur¬ 
ther discussed in Section IV. In that setting, the relevant over¬ 
lap between states |0(G)) and |0'(G')) is entirely encoded 
in their bond matrices G and G' alone. Accordingly, we de¬ 
fine the overlap (0'(G')|0(G)) by comparing G and G' in 
the scalar product that takes into account the gauge difference 
between A and A': 


(0'(G')|0(G)) = 


G 

G'* 


(A12) 


where P is the left dominant eigenvector of the mixed trans¬ 
fer matrix T defined in (All). Similarly, Q is the right 
dominant eigenvector of the mixed transfer matrix built by 
tensors B and B'. P and Q are normalized such that 
(0'(A')|0(A)) = 1. 
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